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We have reconstructed the quantum state of optical 
pulses containing single photons using the method of phase- 
randomized pulsed optical homodyne tomography. The 
single-photon Fock state jl) was prepared using conditional 
measurements on photon pairs born in the process of para- 
metric down-conversion. A probability distribution of the 
phase-averaged electric field amplitudes with a strongly non- 
Gaussian shape is obtained with the total detection efficiency 
of (55 ± 1)%. The angle-averaged Wigner function recon- 
structed from this distribution shows a strong dip reaching 
classically impossible negative values around the origin of the 
phase space. 



a. Introduction States of quantum systems can be 
completely described by their Wigner functions (WF), 
the analogues of the classical phase-space probability dis- 
tributions. Generation of various quantum states and 
measurements of their WFs is a central goal of many ex- 
periments in quantum optics [1-3]. Of particular inter- 
est are quantum states whose Wigner function takes on 
negative values in parts of the phase space. This clas- 
sically impossible phenomenon is a signature of highly 
non-classical character of a quantum state. 

Quantum states containing a definite number of energy 
quanta (Fock states |n)) are paradigmatic in this respect. 
Their WFs exhibit strong negativities and their marginal 
distributions are strongly non-Gaussian (Fig. 1). This 
property reflects the fundamentally non-classical nature 
of these states as carriers of the particle aspect of light. 




FIG. 1. Theoretical phase space quasiprobability density 
(Wigner function) of the single-photon 

state W(X,P) = |(4(X2 + P^) - l) e-2(^'+^'>. 

X — [a + a^)/\/2 and P — [a — a^)/\/2j are normalized 
non-commuting electric field quadrature observables. Sin- 
gle-quadrature probability densities (marginal distributions) 
are also displayed. 

Generation and complete measurements of the Fock 



states' WFs were performed on vibrational states of a 
trapped Be^ ion [y. In the electromagnetic domain, 
Nogues et al. recently reported the measurement of 
the WF of a single-photon state in a superconducting 
microwave cavity at a single point (origin) of the phase 
space. A full characterization of a Fock state of the elec- 
tromagnetic field has not been achieved so far. 

In this paper we present a measurement of the com- 
plete (phase-averaged) Wigner function of the propagat- 
ing single-photon state 1 1) in the optical domain. We per- 
form a direct measurement of the dynamical variables of 
the electromagnetic field, the electric field quadratures, 
whereby their probability distributions are obtained. The 
Wigner function is then reconstructed from the measured 
distributions. This method, homodyne tomography, has 
been established as a reliable technique of reconstructing 
quantum states in the optical domain. Previously, it has 
been applied to classical and weakly non-classical states 
of the light field, such as vacuum, coherent, thermal and 
squeezed states, in the continuous-wave as well as in the 
pulsed regime . 

The main challenge associated with a tomographic 
characterization of the single-photon state is the prepara- 
tion of this state in a well-defined spatio-temporal mode. 
We solve this task by employing conditional state prepa- 
ration on a photon pair born in the process of parametric 
down-conversion The two generated photons are 

separated into two emission channels according to their 
propagation direction (Fig. 2). A single-photon counter 
is placed into one of the emission channels (labeled trig- 
ger) to detect photon pair creation events and to trigger 
the readout of a homodyne detector placed in the other 
(signal) channel Q. 
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FIG. 2. Simplified scheme of the experimental setup 

b. Theory The process of pulsed 2-photon down- 
conversion produces strongly correlated photon pairs. 
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The generated biphoton state can be written as 
( |0, 0} + / dks dkt ^ks,kt) 



, (1) 



where TV is a normahzation constant, k^ and kt denote 
the signal and trigger beam wave vectors, respectively, 
and the function (f>(/cs,/ct) carries the information about 
the amplitude as well as the transverse and longitudinal 
structure of the photon pair generated 

A single-photon Fock state is prepared from \'^) by 
projecting this state onto a photon count event in the 
trigger beam path: 



(2) 



where pt denotes the state ensemble selected by the trig- 
ger and the trace is taken over the trigger states. 

The trigger state ensemble pt is determined by the spa- 
tial and spectral filtering in the trigger channel: 



Pt 



(3) 



where T[kt) is the spatiotemporal transmission function 
of the filter. Note that although \^) represents a pure 
state, Pt and hence ps are statistical mixtures. How- 
ever, if sufficiently tight filtering is applied to the trigger 
channel (so that T(kt) is much narrower than the spatial 
and spectral width of the pump beam) , ps will approach 
a pure single-photon state Q. In this case, the signal 
photons are prepared in a relatively well-defined optical 
mode suitable for homo dyne detection. 

It is important to understand that the "signal beam" 
as shown in Fig. 2 is not an optical beam in the traditional 
sense. The down-converted photons are in fact emitted 
randomly over a wide solid angle. The optical mode of 
the signal state is created non-locally only when a photon 
of a pair hits the trigger detector and is registered. The 
coherence properties of this mode are determined by the 
optical mode of the pump and the spatial and spectral 
filtering in the trigger channel. 

Once the approximation ps of the Fock state is pre- 
pared, it is subjected to balanced homodyne detection. 
The signal wave is overlapped on a beamsplitter with a 
relatively strong local oscillator (LO) wave in the match- 
ing optical mode. The two fields emerging from the 
beamsplitter are incident on two high-efficiency photo- 
diodes whose output photocurrents are subtracted. The 
photocurrent difference is proportional to the value of the 
electric field operator E{9) in the signal mode, 9 being 
the relative optical phase of the signal and LO. 

For each phase 9 one measures a large number N of 
samples of E{9) oc Xe = X cos9 + Psm9, so that their 
histogram (i. e. the marginal distribution) pT{Xg) can be 
determined. The latter is related to the WF as follows: 

pr(Xe) - {Xe\pmce.\Xe) (4) 

/OO 
W{Xcos9- P sin 0,Xsm9 + P cos 9) dP, 
-oo 



where Pmcas is the density matrix of the state being mea- 
sured. The marginal distribution pr(Xe) can be envi- 
sioned as a density projection of the WF W{X, P) onto 
a vertical plane oriented at an angle 9 with respect to the 
plane P — (Fig. 1). From the set of marginal distribu- 
tions pr(Xe) for a large number of phase angles 9 the WF 
of Pmcas can be reconstructed via a procedure similar to 
the one used in medical computer tomography 

In a perfect experiment, Pmeas = |1)(1|7 where the 
single-photon state is in the optical mode which matches 
that of the local oscillator. In reality, various imperfec- 
tions (such as optical losses in the signal arm, inefficient 
photodiodes, dark counts, non-ideal matching between 
the signal and the LO optical modes ||l^) cause an ad- 
mixture of the vacuum |0) to the measured state, so that 



Pn 



77|l)(l| + (l-,/)|0)(0|, 



(5) 



rj being the measurement efficiency. It is remarkable 
that all these effects act upon Pmcas in a similar way, 
so that their effect can be expressed in a single number 
■q which is a product of efficiencies associated with indi- 
vidual parts of the setup. The value of 77 is crucial for 
this experiment as it strongly influences the shape of the 
measured marginal distributions and the reconstructed 
Wigner function (Fig. 3) 
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FIG. 3. Effect of the non-perfect measurement efficiency 
rj on the marginal distribution (a) and the reconstructed WF 
(b). For the WF, cross-sections by the plane P = Q are shown. 
Negative values require rj > 0.5. 

In our experiment we used a simplified scheme in 
which the phase 9 varied randomly, so that we only 
measured a single phase-randomized marginal distribu- 
tion prav(-'f) = (pr(Xe))e. This does not change the 
measurement result for quantum states with rotationally 
symmetric Wigner functions such as those described by 
Eq. d). The phase -averaged WF W{R) is obtained from 
prav(^) via the Abel transformation |^,|ll[: 



W{R) 



dprav(X) 
dX 



- R^)-^/^dX. (6) 



From the phase-randomized marginal distributions one 
can also directly infer diagonal elements p„„ of the state 
density matrix in the Fock basis. 



(7) 



Pnn = TI" / prav(Ar) /mi(X) dX , 



where fnn{X) are the amplitude pattern functions [T^ ] 
which are independent of the optical state being sampled. 
The statistical uncertainty of the reconstructed pnn is 



2 



i'^L) =Nj_ P^av(^) {^fnn{X)f dX , (8) 

where N is the total number of field samples acquired. 

c. Experimental Setup We employed a mode-locked 
Ti:sapphire laser (Spectra Physics Tsunami) in combi- 
nation with a pulse picker to obtain transform-limited 
pulses at 790 nm with a repetition rate 816 kHz and a 
pulse width of 1.6 ps. Most of the radiation was fre- 
quency doubled in a single pass through a 3-mm LBO 
crystal yielding 100 /iW at 395 nm and then passed on 
to a 3-mm BBO crystal for down-conversion. 

Down-conversion occurred in a type-I frequency- 
degenerate, but spatially non-degenerate configuration, 
with 790-nm photon pairs emitted at angles ±6.8° with 
respect to the pump beam. The BBO crystal was cut at 
6 = 35.7°, (p — 0°, so that the direction of the walk-off 
of the 395 nm pump beam inside the crystal coincided 
with the direction of the signal beam so as to minimize 
distortions of the signal spatial mode ( "hot spot" down- 
conversion, ) . The short crystal lengths allowed us to 
avoid group-velocity mismatch effects which would have 
complicated temporal mode matching to the LO pulse. 

The trigger photons passed through a spatial filter and 
a 0.3-nm interference filter centered at the laser wave- 
length. They were then detected by an EG&G SPCM- 
AQ-131 single-photon detector (quantum efficiency 60%, 
dark count rate < 15 s~^) at a rate of about 0.25 s~^. 
Such a low pair production rate made the effect of Fock 
states with n > 1 negligible. Precise (within 0.6 ns) gat- 
ing of the count events with the laser pulses allowed us 
to eliminate most of the dark counts, thereby reducing 
their contribution to about 2% of all trigger events. 

We used a small fraction of the original optical pulses 
from the pulse picker — split off before the frequcncy- 
doubler — as the local oscillator for the homodyne sys- 
tem. Achieving a good spatial and temporal mode match- 
ing between the LO and the photons in the signal channel 
constituted a major challenge in this experiment due to 
extremely low intensity of the field in the signal mode. 
To this end, a fraction of the laser output power was di- 
rected into the BBO crystal from the back along the trig- 
ger beam path so that it passed through the spatial filter 
in the trigger channel. Inside the crystal these alignment 
pulses were temporally and spatially overlapped with the 
pump to produce a difference frequency (DFG) emission 
into an optical mode which modeled, to a good precision, 
that of the conditionally prepared signal photons. This 
mode was then matched to that of the local oscillator by 
observing an interference pattern between the two beams. 
A visibility on the level of w = 83 ± 1% was reached. 

The method of conditional state preparation also es- 
tablished special requirements for the homodyne detector 
electronics. The detector needed to resolve quantum shot 
noises of individual laser pulses at a 0.8-MHz repetition 
rate. Details on design and performance of the homodyne 
system developed will be published elsewhere Q . 
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FIG. 4. Experimental results; a) raw quantum noise data 
for the vacuum (left) and Fock (right) states along with their 
histograms corresponding to the phase-randomized marginal 
distributions; b) diagonal elements of the density matrix of 
the state measured, c) reconstructed WF which is negative 
near the origin point. The measurement efficiency is 55%. 

d. Results and discussion In a 14-hour experimental 
run about 200,000 vacuum state and 12,000 Fock state 
samples were acquired (Fig. 4a). Both data sets were 
then binned up to obtain their statistical distributions. 
A Gaussian distribution was fit to the vacuum state noise 
spectrum by varying its X-scale and point of origin. 

The best fit parameters of the vacuum state were used 
to scale the Fock state data. The latter was then fit 
by the theoretical marginal distribution of the ensemble 
(^) to find the measurement efficiency rj. The best fit 
efficiency value was rj = 0.55. Using the Abel transform 
(^ the phase-randomized WF of the observed quantum 
state was reconstructed (Fig. 4c). As expected, it exhibits 
negativity around the origin point, with a minimum value 
IF(0,0) = -0.062. 

The diagonal elements of the density matrix have been 
evaluated along with their statistical errors by apply- 
ing the quantum state sampling method as defined by 
Eqs. (|^, ^) directly to the rescaled raw data. We found 
pii — 0.553 ± 0.013 in agreement with the value of rj 
obtained by fitting the marginal distribution (Fig. 4b). 
Substituting this quantity into Eq. (||) and calculating 
the corresponding value of the WF we find W{0, 0) = 
—0.067 ±0.016, in agreement with the above value deter- 
mined from the marginal distribution. The uncertainty 
of ±0.016 gives an estimation for the accuracy to which 
the negative value of the Wigner function has been de- 
termined in this experiment. 

The values of poo and pn for the vacuum state were 
found to be 0.9975 ± 0.0029 and 0.0021 ± 0.0032, respec- 
tively. These quantities being equal to their ideal values 
within the statistical errors indicates that no unknown 
conditions (e.g. mechanical vibrations) are present that 
might bring the apparent measurement efficiency above 
its actual value. The systematic experimental errors (e.g. 
technical noise in the homodyne detector) were insignif- 
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icant and would result in decrease of the measured effi- 
ciency. It is thus impossible that neglecting these errors 
could bring about a too optimistic estimate of r/. 

What are the main factors reducing the measurement 
efficiency? The non-perfect match between the optical 
modes of LO and DFG fields diminishes the efficiency 
by a factor of = 0.69 ± 0.02. The fact that the DFG 
wave does not perfectly mimic the conditionally prepared 
mode of the single photon causes an additional reduction 
by xO.95. A factor of 0.90 arises from losses in the sig- 
nal beam path and non-perfect quantum efficiency of the 
photodiodes in the homodyne detector. Finally, a 2% 
reduction occurs due to false trigger count events. Com- 
bining all the above factors we obtain the upper limit 
estimate of the quantum efficiency as (57 ±2)%, which is 
in good agreement with the experimental value of 55.3%. 

e. Conclusion and outlook We have reconstructed the 
phase-averaged Wigner function and the density matrix 
diagonal elements of an optical single-photon Fock state 
1) with a total measurement efficiency of 55.3 ± 1.3% 
using the method of phase-randomized pulsed optical ho- 
modyne tomography. The reconstructed WF is of non- 
Gaussian shape and exhibits negative values around the 
origin of phase space, reflecting the strongly non-classical 
character of the state |1) as a particle state of the light 
field. Single-photon Fock states were prepared in a well- 
defined electromagnetic mode by conditional measure- 
ments on photon pairs created in the process of para- 
metric fluorescence. The measurement technique and er- 
ror analysis were checked by performing a simultaneous 
measurement on the vacuum state. Major experimental 
inefficiency factors have been identified and quantified. 
This experiment represents the first quantum tomogra- 
phy measurement of a highly nonclassical state of the 
electromagnetic field. 

Relatively straightforward modifications of the setup 
would allow tomography measurements of displaced Fock 
states and 1-photon added coherent states Of 
special interest is the entangled state |0, 1) — |1, 0) gener- 
ated when a single photon is incident on a beamsplitter. 
This state can also be used to demonstrate nonlocality of 
a single photon In a more distant future, arbitrary 
quantum states of the light field might be generated via 
repeated down-conversion p^ ]. 
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